Renormalization of Molecular Electronic Levels at Metal-Molecule Interfaces 
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The electronic structure of benzene on graphite (0001) is computed using the GW approximation 
for the electron self-energy. The benzene quasiparticle energy gap is predicted to be 7.2 eV on 
graphite, substantially reduced from its calculated gas-phase value of 10.5 eV. This decrease is caused 
by a change in electronic correlation energy, an effect completely absent from the corresponding 
Kohn-Sham gap. For weakly-coupled molecules, this correlation energy change is seen to be well 
described by a surface polarization effect. A classical image potential model illustrates trends for 
other conjugated molecules on graphite. 
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There is renewed interest in using organic molecules as 
components in nanoscale electronic and optoelectronic 
devices 0, 0] , and thus a critical need has emerged for 
improved knowledge and control of charge transport phe- 
nomena in organic molecular assemblies |3| . Understand- 
ing transport across the interface between the active or- 
ganic layer and the metallic electrode has proved partic- 
ularly challenging, especially in the single-molecule limit. 
Fundamentally, charge transport is controlled in such 
systems by the electronic coupling of frontier molecu- 
lar orbitals to extended states in the electrode, and the 
energetic position of these orbitals relative to the con- 
tact Fermi level. Several recent measurements of or- 
ganic thin films, self-assembled monolayers (SAMs), and 
single- molecule junctions have emphasized the important 
role of Coulomb interactions between the added hole or 
electron in the frontier orbitals and the metal substrate 
milSSIilillSIIS- However, most theoretical cal- 
culations of transport through organic molecules have 
continued to rely on some implementation of density 
functional theory (DFT) or semiempirical one-particle 
Hamiltonians Q . The limitations of DFT for describing 
excited-state energies are well known |l2j . and implica- 
tions for a DFT-based theory for nanoscale conductance 
have been recently discussed 0|. 

When a molecule is brought in contact with a metal, 
several physical effects will influence its ionization level 
(highest occupied molecular orbital, HOMO) and affin- 
ity level (lowest unoccupied molecular orbital, LUMO). 
First, the self-consistent interaction between molecule 
and surface will rearrange the electron density and mod- 
ify the alignment of frontier orbital energies. Second, 
electronic coupling to extended states in the metal will 
further shift orbital energies and broaden discrete molec- 
ular levels into resonances. Finally, the Coulomb interac- 
tion between the added hole or electron associated with 
the ionization or affinity level will result in a polariza- 



tion of the metal substrate. This additional correlation 
energy further stabilizes the added hole or electron, re- 
ducing the gap between affinity and ionization levels as 
illustrated in Fig. 1. An accurate DFT-based approach 
should correctly capture the first effect j although the 
use of DFT to calculate the width of resonances is under 
debate 13] . Importantly however, the surface polariza- 
tion response, as we show here, is completely absent from 
frontier orbital energies computed in DFT. 
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FIG. 1: Schematic energy level diagram indicating polariza- 
tion shifts in the frontier energy levels (ionization and affinity) 
of a molecule upon adsorption on a metal surface. 

In this Letter, we compute the electronic excited states 
for a clean, weakly-coupled system consisting of an aro- 
matic molecule (benzene) physisorbed on the graphite 
(0001) surface. Electronic correlations are included di- 
rectly within a first-principles many-electron Green func- 
tion approach The electron self-energy is calcu- 
lated from first principles within the GW approximation 
(GWA) [l^ using a methodology 0| that has proved 
accurate for a wide range of systems |0. While more 
generally including dynamical electronic correlation, the 
GWA is well known to include static, long-ran ge i mage 
potential effects for an electron near an interface • Us- 
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ing this theoretical approach, we predict a strong renor- 
mahzation of the electronic gap of the benzene system 
(relative to its molecular gas-phase value) when it is ph- 
ysisorbed on a graphite (0001) surface. The change in 
the electron correlation energy on adsorption can be un- 
derstood as a polarization effect in this case. An image 
potential model is used to illustrate trends for other aro- 
matic molecules. 

Equilibrium geometries of molecular benzene in the 
gas-phase, condensed in a bulk crystalline phase, and ph- 
ysisorbed on graphite (0001) are determined using DFT 
within the local density approximation (LDA). Norm- 
conserving pseudopotentials |20| are used with a plane- 
wave basis for the electron wavefunctions (80 Ry cutoff) 
for structural relaxations. The surface is modeled with 
a 3x3 supercell containing 4 layers of graphite, a single 
benzene molecule, and the equivalent of 7 layers of vac- 
uum. The theoretical in-plane bulk lattice parameter is 
used (a=2.45 A, c=6.62 A). In the most stable site for 
adsorption, benzene rests flat on the surface centered on 
a three-fold site 3.25 A above a substrate carbon atom, 
in agreement with a previous study [2]|. For compar- 
ison, benzene is also considered in an upright position, 
centered above a hollow site with its closest hydrogen 
atom 2.21 A from the surface. Solid crystalline benzene 
has an orthorhombic unit cell containing four molecules 
(Pbca); the atomic positions within the unit cell are op- 
timized keeping the lattice parameters a, b and c fixed 
to their experimental values of of 7.44, 9.55 and 6.92 A 
respectively ^ ■ The gas phase is modeled using a cubic 
supercell (a=13.22 A). For each system, matrix elements 
of the self-energy operator are evaluated using a 50 Ry 
energy cutoff for the electronic wavefunctions, a 6 Ry 
cutoff for the momentum-space dielectric matrix, and a 
2.9 Ry cutoff for the sum on the virtual states. This 
choice of parameters results in quasiparticle energy gaps 
converged within ~ 0.2 eV. 

The electron addition and removal energies of a ben- 
zene molecule in the gas-phase, calculated in the present 
GW approach, result in a HOMO-LUMO (quasiparticle) 
gap of 10.51 eV. This value agrees well with an inde- 
pendent GW calculation 1231 . total energy difference cal- 
culations based on DFT |24L Eij i. and experiment [2fij |. 
By contrast, the Kohn-Sham gap (within LDA) is 5.1 
eV, substantially smaller. The electronic structure of 
benzene on the graphite (0001) surface along the F-K' 
direction is shown in Fig. 2. Comparing the surface- 
projected band structures (shaded regions) in Fig. 2(a) 
and 2(b), the quasiparticle bandwidth increases by about 
15% relative to LDA, in agreement with previous works 
I27I . The bold horizontal lines interpolate between 
the benzene HOMO and LUMO states computed at F 
and K'; the filled circles at these high-symmetry points 
indicate states with significant weight on the molecule. 
For physisorbed benzene, the Kohn-Sham (LDA) gap is 
5.05 eV throughout the zone, unchanged from the cor- 
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FIG. 2: Calculated frontier orbital energy levels (heavy blue 
lines) with the indicated energy gap (red arrows) for benzene 
adsorbed flat on the graphite surface, plotted against the pro- 
jected surface band structure of the graphite: (a) DFT (LDA) 
energies, (b) GW quasiparticle energies, (c) GW quasiparti- 
cle energies of benzene in the gas phase. Inset in (a) shows a 
model of the adsorption geometry. 



responding LDA gas-phase value. Relative to the LDA 
value, the quasiparticle gap of the molecule flat on the 
graphite surface is much larger, 7.35 eV. However, the 
predicted quasiparticle gap is substantially smaller than 
the gas-phase value of 10.51 eV. 

Table 1 summarizes the calculated HOMO-LUMO 
gaps of benzene in the four environments considered in 
this study. Interestingly, the LDA gaps are identical for 
all environments. In contrast, the GW self-energy cor- 
rections exhibit noticeable variation. To understand this 
variation, we analyze the self energy change relative to 
the gas-phase. The change AE for each frontier level 
is decomposed into Coulomb- hole (AEch), screened- 
exchange (AEsx), and bare exchange or Fock (AEx) 
contributions. We find that ASch is nearly equal for 
the occupied and empty frontier states, and that ASx is 
quite small (0.1-0.2 eV). Interestingly, the screened ex- 
change term that is responsible for most of the difference: 
for the HOMO, we observe ASsx ~ -2AI;ch, while for 
the LUMO AEsx ^ 0. Put together, the change in cor- 
relation energy (AScorr — ASch + ASsx) reported in 
Table 1 turns out to be nearly symmetric between the 
ionization and affinity levels for the benzene in each en- 
vironment studied. This result is qualitatively the same 
as that obtained from the derivation of the image poten- 
tial effect for an electron near a metal surface [l9j. 

To develop a more detailed model of our self energy re- 
sults, we recognize that the benzene frontier orbitals are 
only weakly coupled to the environment. When the or- 
bitals have negligible overlap with the metal, the correc- 
tion to the molecular self-energy operator upon adsorp- 
tion depends only on the change in the screened Coulomb 
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TABLE I: Benzene HOMO-LUMO gaps in the gas phase, 
crystal phase, and adsorbed on the graphite surface (fiat and 
perpendicular) . First and second lines are Kohn-Sham (LDA) 
and quasiparticle (GW) gaps. (For the crystal, we average 
over the vr and tt* manifolds.) Third and fourth lines are 
calculated changes in correlation energy for the HOMO and 
LUMO, relative to the gas phase, determined from the full 
GW calculations and from an image potential model. Ener- 
gies are in eV. 



Gas Flat Perp Crystal 
phase graphite graphite phase 
AEgap (LDA) 5l6 KoE 5T1 K07 
ASgap (GW) 10.51 7.35 8.10 7.91 
AEcorr 1.45, -1.51 1.18, -1.17 1.16, -1.15 
AScorr (Model) 1.50, -1.43 0.97, -0.96 



interaction W, i.e. 

occ 

AEsx(r,r';i?) = ^0j(r)^*(r')AW^(r,r';£;-i?,), (1) 
j 

where are molecular wavef unctions and Ej their 
eigenvalues. A corresponding expression exists for the 
Coulomb-hole term. For sufficiently large metal-molecule 
separations, AW is smooth and slowly-varying over the 
spatial extent of the molecular orbitals, and only the self- 
term contributes to single-particle matrix elements of Eq. 
(1). Then the change in correlation energy from the sur- 
face can be reduced to 

Ai?HOMo — (^HOMolASsx + ASchI^homo) 
^ 2Phomo — ^HOMO = -Phomo 

and 

ASlumo ~ (0LUMo|AI]sx + AEchI^lumo) 

^ — PlUMO — —^LUMOi 

where P is the static polarization integral 

P^^-\ J J drdrV, (r)0*(r')AW^(r,r')0,(r')0*(r). 

(4) 

For benzene on graphite, the full GW calculations indi- 
cate that dynamical effects make a negligible contribu- 
tion to AScorrj and that the self-term accounts for more 
than 90% of AEcorr, supporting the simplified picture of 
Eqs. (2-4). 

Significant further simplification is achieved if an im- 
age potential model is sufficient for APF(r,r). In Fig. 3, 
the screening potential through the molecular adsor- 
bate is illustrated by considering AH^(r, r) = AT4cr(r, r), 
where Vl5cr(r, r) results from the screening response to 
the added electron (or hole) The difference be- 

tween the adsorbed molecule and isolated molecule, 
AT4cr(r,r), is compared with the image potential model. 
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FIG. 3: Static screening potential, I4cr(r, r), plotted along a 
line through benzene adsorbed flat on graphite. Thin, solid 
(black) curve is the total Vscr(r,r) for the metal- molecule - 
system; the thick, solid (red) curve is A14cr(r,r), the change 
upon adsorption; the heavy, dashed (blue) curve is the image 
potential model relative to the image plane (light, vertical 
dashed line). Inset: physical model, scaled to the axis of 
the plot, including an isosurface plot of a frontier benzene tt 
orbital. 

AVscrlr, r)-^l/4|z — zo\. The image plane position zq is 
explicitly determined, in a separate calculation |2^, to 
be 1 A beyond the last surface plane for our graphite 
slab. From Fig. 3 it can be seen that, over the spa- 
tial range of the molecular orbital, an image potential 
captures the main effect. Using the value of the image 
plane position computed above and the frontier orbitals, 
-Phomo and Plumo for benzene flat and perpendicular 
on the graphite surface are calculated using the image 
potential model. As shown in Table 1, the image poten- 
tial model is quite accurate for the flat case and captures 
most of the effect for the perpendicular case (within 0.2 
eV). The simple image potential model neglects the in- 
ternal screening response of the molecule to the polariza- 
tion of the metal surface. While small for a flat molecule 
oriented parallel to a surface, a significant molecular po- 
larizability perpendicular to the metal surface leads to 
an additional energy gain, increasing the Pj. 

Provided that molecular resonances are well separated 
from the metal Fermi energy, the polarization model for 
AE should be broadly applicable. To illustrate the im- 
pact of the change in correlation energy, we use the image 
model to predict the renormalized gaps for members of 
the acene series and coronene adsorbed flat on graphite 
(Table 2). For the larger molecules in the series, the 
change in gap becomes dramatic, e.g. the pcntacene gap 
is predicted to diminish by nearly a factor of two on a 
graphite surface. 

The role of geometry and morphology on changes in po- 
larization energy in organic systems can be subtle |30j. 
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TABLE IL For selected molecules, measured gas-phase ioniziation energies and electron affinities are combined with an 
image model for polarization energies to predict adsorbate HOMO-LUMO gaps for molecules flat on graphite (all in eV). 



Molecule 


Expt IP 


Expt EA 


Gas-phase gap 


ftoMO 


PhVMO 


Adsorbate gap 


Naphthalene 


8.14 


-0.20 


8.34 


1.41 


1.39 


5.54 


Anthracene 


7.44 


0.53 


6.91 


1.32 


1.30 


4.29 


Tetracene 


6.97 


0.88 


6.09 


1.24 


1.23 


3.62 


Pentacene 


6.63 


1.39 


5.24 


1.18 


1.18 


2.88 


Coronene 


7.29 


0.47 


6.82 


1.19 


1.17 


4.46 



but the impact has been directly measured for organic 
films on metal substrates using photoemission and in- 
verse photoemission 0. Adsorbate frontier orbital en- 
ergies can be probed directly by STM, provided the 
HOMO-LUMO gap is small enough for the resonant tun- 
neling regime to be experimentally accessible 0, l3lj . 
From Table 2, tetracene and pentacene are within typi- 
cal measurement range (±2.5 V), while anthracene and 
coronene are marginal. In a recent study of pentacene ad- 
sorbed on ultrathin NaCl on Cu(lll) jllj, gaps of 3.3, 4.1 
and 4.4 eV are observed for NaCl thicknesses of one, two 
and three monolayers, respectively. Our modeled value 
of 2.9 eV for direct adsorption on the graphite surface 
fits well with this progression. 

In conclusion, we find that the correlation contribution 
to the frontier molecular orbital energies depends sensi- 
tively on environment. In the examples studied here, 
the change in correlation energy is dominated by a po- 
larization effect. The impact of electrode surface polar- 
ization on spectroscopic measurements must be carefully 
assessed for each metal-molecule system. For organic 
films or SAMs on a metal, the polarization contribution 
from neighboring molecules can also be quite significant. 
For molecular systems where the frontier orbitals have 
stronger electronic coupling to the metal and the resul- 
tant resonances overlap with the Fermi energy, the role 
of dynamical charge transfer is expected to be consider- 
able. Future investigations must address the nature of 
additional contributions to the self energy in the event of 
stronger coupling. 
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